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SUMMARY 


+*-  this  paper  we-  study-  the  probabilistic  behavior  of  a two-echelon 

inventory  system  consisting  of  a depot  and  a set  of  bases.  Primary 

demands  occur  at  the  bases  for  a single  unit  at  a time.  Whenever  a 

base's  inventory  position  reaches  the  reorder  point  (s),  the  base 

orders  sufficient  inventory  from  the  depot  to  raise  the  base's  inventory 

position  to  S.  Similarly,  the  depot  places  an  order  to  its  supplier 

when  its  inventory  position  reaches  its  reorder  point;  the  order  is  for 

the  number  of  units  required  to  raise  the  depot's  inventory  position 

to  a prespecified  level.  Thus  all  locations  follow  a continuous  review 

(S,s)  policy.  All  excess  demand  is  assumed  to  be  backordered. 

Tlie. 

Our  main  objective  is  to  derive  the  probability  distribution  for 
the  number  of  backordered  units  at  a base  at  an  arbitrary  point  in  time 
given  an  item  follows  a known  (S,s)  policy  at  each  base  and  the  depot. 
The  demand  process  at  each  base  is  assumed  to  be  a stationary  Poisson 

process.  The  analysis  is  carried  out  for  two  cases.  In  the  first  case, 

it  ii  . 

we-  assumed  the  system  consists  of  a large  number  of  bases;  in  the  second 

){  i'l  . . 

case,  we  assumed  there  are  two  bases  in  the  system.  To  simplify  the  dis- 
cussion, ■we  assume^ in  both  cases  that  all  bases  follow  the  same  (S,s) 


policy . 


I.  INTRODUCTION 


In  this  paper  we  study  the  probabilistic  behavior  of  an  inventory 
system  consisting  of  a depot  (the  first  echelon),  and  a set  of  bases 
(the  second  echelon).  Primary  demands  occur  only  at  the  bases.  Each 
demand  at  a base  is  for  a single  unit  of  an  item.  When  a base's  inven- 
tory position — on-hand  plus  on-order  minus  backorders — reaches  the 
reorder  point  (s),  the  base  orders  enough  inventory  from  the  depot  to 
raise  the  base's  inventory  position  to  a pre-established  level  (S). 
Similarly,  whenever  the  depot  inventory  position  reaches  its  reorder 
point,  the  depot  orders  a sufficient  quantity  from  its  supplier  to 
raise  its  inventory  position  to  a given  value.  Thus  each  location  in 
the  system  follows  a continuous  review  (S,s)  policy.  Furthermore, 
demands  occurring  when  a base  is  out  of  stock  are  assumed  to  be  back- 
ordered . 

Numerous  real  inventory  systems  operate  in  this  manner.  For  example, 
the  Air  Force  maintains  bases  throughout  the  world.  Demand  for  spare 
parts  occurs  at  these  bases;  the  bases  in  turn  are  resupplied  for  most 
components  by  an  Air  Force  or  other  Department  of  Defense  depot.  For  each 
component  there  is  only  one  depot  that  resupplies  a base.  Each  location 
in  this  system  follows  a continuous  review  (S,s)  type  inventory  policy. 

Our  main  objective  is  to  develop  the  probability  distribution  for 
the  number  of  backordered  units  at  a base  at  an  arbitrary  point  in  time 
for  a specific  item  given  a particular  (S,s)  policy  for  each  base  and 
the  depot.  The  approach  we  take  can  also  be  used  to  determine  the 
probability  distribution  for  on-hand  inventory  at  a base  at  an  arbitrary 


point  in  time.  We  assume  the  demand  process  at  each  base  is  a stationary 
Poisson  process.  The  analysis  is  carried  out  for  two  cases. 

In  the  first  case,  we  assume  the  item  is  used  at  a large  number  of 
bases.  For  example.  Xerox  has  many  spare  parts  stocked  by  most  of  the 
approximately  8000  people  who  repair  Xerox  machines.  Each  person,  cor- 
responding to  a base,  is  normally  resupplied  by  one  Branch  warehouse, 
corresponding  to  a depot.  Roughly  90  people  are  resupplied  by  a single 
Branch  warehouse.  To  simplify  the  discussion,  we  will  assume  the  demand 
process  is  the  same  at  each  base.  Furthermore,  the  costs — carrying  costs, 
shortage  costs,  and  fixed  and  variable  ordering  costs--and  depot-to-base 
transportation  times  are  assumed  to  be  the  same  for  each  base.  If  each 
base  follows  an  (S,s)  policy,  they  should  all  select  the  same  values 
for  S and  s;  we  assume  that  they  do.  (The  analysis  we  present  can 
be  modified  using  methods  similar  to  those  used  in  Chapter  5 of  reference  8 
to  study  the  case  where  demand  distributions,  costs,  or  transportation 
times  are  not  the  same  at  all  bases). 

Next,  we  examine  a system  consisting  of  two  identical  bases  and  a 
depot.  The  demand  processes,  costs,  and  transportation  times  are  the  same 
at  both  bases.  Each  base  is  assumed  to  follow  the  same  (S,s)  policy. 
Again,  we  could  easily  extend  the  results  to  the  case  where  there  are 
more  than  two  bases;  however,  for  ease  of  discussion  we  will  describe 
in  detail  only  the  two  base  case.  As  an  example  of  this  case,  the 
Air  Force  operates  FB-111  aircraft  mainly  at  two  locations  at  which 
approximately  the  same  removal  rates  are  observed  for  most  components. 

An  excellent  survey  of  many  papers  on  multi-echelon  inventory  systems 
is  given  by  Clark  [1].  Most  of  these  papers  deal  with  the  periodic 
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review  situation.  Our  work  specifically  extends  the  results  given  by 
Simon  L6]  and  Kruse  and  Kaplan  [4].  Their  papers  pertain  to  two-echelon 
systems  in  which  the  depot  follows  a continuous  review  (S,s)  policy, 
but  the  bases  are  restricted  to  a continuous  review  (S,S-1)  policy. 

Before  proceeding  with  the  analysis,  we  pause  to  comment  on  the 
motivation  for  this  research.  This  paper  is  part  of  a larger  study 
which  examines  the  impact  on  cost  and  performance  of  following  different 
inventory  policies  in  several  real  multi-item,  multi-echelon  systems  [6]. 
The  results  developed  here  will  be  used  to  evaluate  the  inventory  perfor- 
mance obtained  using  different  values  of  the  policy  variables  for  each 
item  at  each  location  in  one  multi-item  system.  In  particular,  the 
probability  distribution  derived  in  this  paper  will  be  used  to  find  the 
optimal  values  for  the  policy  variable  for  each  item  in  a real  system 
when  the  objective  is  to  select  the  inventory  levels  that  minimize  total 
expected  base  backorders  subject  to  a constraint  on  system  inventory 
investment.  These  values  will  then  be  compared  to  those  obtained  using 
several  approximation  methods.  General  conclusions  will  be  drawn  in 
that  study  concerning  the  quality  and  appropriateness  of  the  several 
approximation  methods  for  determining  policy  variable  values  for  items 
having  different  demand  and  cost  characteristics  when  different  budget 
restrictions  are  imposed  on  the  system. 


4 


II.  SITUATION  ONE-- A LARGE  NUMBER  OF  BASES 

In  this  section  we  develop  the  probability  distribution  for  the 
number  of  units  backordered  at  a particular  base  at  a random  point  in 
time  when  the  number  of  bases  in  the  system  is  large.  The  specific 
assumptions  we  make  are  as  follows. 

1.  All  bases  and  the  depot  follow  a continuous  review  (S,s) 

policy,  and  all  the  bases  use  the  same  values  for  S and  s. 

2.  No  partial  fill  of  base  orders  is  permitted;  that  is,  all 

S-s  units  must  be  shipped  simultaneously  from  the  depot  to  satisfy  a 
base  order. 

3.  A stationary  Poisson  process  generates  demand  at  each  base. 
Furthermore,  the  demand  rate  A is  the  same  at  each  base. 

4.  Demands  occurring  when  a base  has  no  on-hand  inventory  are  back- 
ordered . 

5.  The  depot-to-base  transportation  time  T is  constant  and  the 
same  for  all  bases . 

6.  The  depot  resupply  time  R is  constant. 

7.  The  depot  reorder  point  r^  is  greater  than  or  equal  to  minus 

one;  the  base  reorder  points  are  also  greater  than  or  equal  to  minus  one. 

8.  The  number  of  bases  in  the  system  is  large. 

We  now  discuss  the  implications  of  these  assumptions. 

Since  the  demands  at  the  bases  are  generated  one-at-a-time , the  (S,s) 


policy  is  the  same  as  a (Q,r)  policy,  where  r indicates  the  reorder 
point  and  Q = S - s.  Hence,  we  will  refer  to  the  base  policy  as  a (Q,r) 
policy . 
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Next,  we  make  an  observation  about  the  depot's  inventory  position. 

By  assumption,  a partial  fill  of  a base's  order  is  prohibited.  Furthermore, 
we  observed  that  our  assumptions  imply  that  all  bases  should  follow  an 
identical  (Q,r)  policy.  It  then  follows  that  the  depot  inventory 
position  should  always  be  a multiple  of  Q.  That  is,  there  would  be  no 
advantage  to  having  the  inventory  position  be  other  than  a multiple  of 
Q since  extra  holding  costs  would  be  incurred  without  improving  the 
chance  of  filling  orders.  Thus  the  depot  always  orders  in  multiples  of 
Q;  that  is,  the  depot  reorder  quantity,  Q^,  can  be  expressed  as 

Qd  = M.Q 


for  some  M = 1,2,...  . 

One  can  show  that  Assumptions  1 and  3 imply  that  the  probability  dis- 
tribution for  the  depot  inventory  position  is  uniformly  distributed  over 
J = {rD+Q,  rD+2Q, . . . »rD+QD}  t5J- 

Our  last  observation  is  a consequence  of  Assumption  8.  Since  the 
number  of  bases  is  large  and  the  successive  times  between  placing  orders 
by  a base  for  depot  resupply  form  a renewal  process,  the  order  arrival 
process  at  the  depot  can  be  accurately  approximated  by  a Poisson  process 
[3]. 

Before  proceeding  with  the  derivation  of  the  probability  distribution 
of  base  backorders,  let  us  introduce  the  nomenclature  used  in  the  remainder 


of  the  paper: 
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Ig  represents  the  depot  inventory  position  at  time  t-T-R, 

I1  represents  base  j's  inventory  position  at  time  t-T-R, 

2 

I represents  base  j's  inventory  position  at  time  t-T, 

G represents  the  number  of  orders  placed  by  all  bases  other 

than  base  j for  depot  resupply  during  (t-T-R, t-T], 

D represents  the  number  of  demands  occurring  during  (t-T-R, t-T] 
at  base  j , 

% 

represents  the  number  of  demands  occurring  during  (t-T,t] 
at  base  j , 

V represents  the  number  of  satisfied  orders  placed  by  base  j 
on  the  depot  during  (t-T-R, t-T]  --  the  orders  are  placed 
during  (t-T-R, t-T]  and  received  at  base  j prior  to  t, 

t \ -x  n 

a(x,n)  = e x /n! 5 

B(t)  represents  the  number  of  backorders  existing  at  base  j at 
time  t, 

U represents  the  number  of  orders  placed  on  the  depot  during 
(t-T-R, t-T]  by  base  j that  are  unfilled  at  time  t, 

represents  the  arrival  rate  of  orders  at  the  depot  from  all 
bases  except  base  j measured  in  orders  per  day. 


and  y 
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The  effect  of  time  at  base  j is  displayed  in  Figure  1. 


Inventory 
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Demand 

Inventory 

Position 

Demand 
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1 

1 

1 
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1 

t-T 

| t line 

t 

4 Depot  Resupply  Time ► < — Depot-to-Base  Transportation-* 

Time 

Time  Sequence  of  Events  at  Base  j 
Figure  1 

We  will  consider  the  steady  state  behavior  of  this  system.  In  partic- 
ular, our  interest  is  in  the  steady  state  distribution  of  the  number  of 
backorders  at  a base;  that  is,  we  want  to  find  P(B(t)  = b)  for 
b = 0,1,...  . To  compute  it,  observe  that  any  order  placed  by  base  j 
prior  to  time  t-T-R  has  been  satisfied  by  time  t.  However,  some  orders 
placed  on  the  depot  during  the  interval  (t-T-R, t-T]  by  base  j may  not 
be  satisfied  by  time  t.  This  could  cause  backorders  to  exist  at  time  t 
for  some  of  the  units  demanded  during  (t-T-R, t-T].  Observe  that  the 

inventory  available  to  satisfy  demands  during  the  interval  (*-T,t]  is 
2 

I -U  • Q.  If  the  demand  at  base  j during  (t-T,t]  exceeds  this  amount 
or  this  quantity  is  negative,  then  backorders  will  exist  at  time  t at 
base  j . 

To  compute  P(B(t)  = b)  we  will  first  compute  P(U  = u).  We 
examine  two  cases.  In  the  first  case,  which  we  call  Case  A,  the  total 
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demand  on  the  depot  during  the  interval  (t-T-R,t-T]  does  not  exceed 
the  available  depot  inventory  (IQ),  so  that  U = 0.  In  the  second 
case,  which  we  call  Case  B,  the  total  depot  demand  during  (t-T-R,t-T] 
exceeds  IQ,  so  that  U may  be  positive. 

First  consider  Case  A,  when  the  number  of  units  ordered  from  the 
depot  during  the  interval  (t-T-R,t-T]  is  not  larger  than  the  available 
depot  inventory.  Note  that  [(D-I^+Q+rJ/Q]  represents  the  number  of 
orders  placed  by  base  j during  the  interval  (t-T-R,t-T],  where  [k] 
represents  the  greatest  integer  less  than  cr  equal  to  k.  Thus, in  this 
case,  G + [ ( D-I *>Q+r )/Q]  < Iq/Q.  As  we  observed,  in  this  situation  all 
orders  placed  by  base  j on  the  depot  prior  to  time  t-T  will  be 
satisfied  by  time  t,  and,  therefore,  P(U  = o|Case  A)  = 1.  Thus  all 
shortages  existing  at  time  t at  base  j are  due  to  demands  placed  at 
base  j during  (t-T,t], 

Then 

r+Q  ^ 

P(B(t)  = b;  Case  A)  = £ P(D  = i+b|Case  A;  I = i)  • P(I  = i;  Case  A) 

i=r+l 

r£Q  >v  2 

= l P(D  = i+b)  • P(I  = i;  Case  A) 

i=r+l 

r+Q  9 

= l a (XT,  i + b)  • p(r  = i;  Case  A),  for  b >_  1,  (Eq.  1) 

i=r+l 


since  D does  not  depend  on  events  prior  to  time  t-T.  We  will  next 

find  P(I2  = i;  Case  A). 

2 1 

Suppose  I = i,  I = k,  and  k < i.  Then  demand  at  base  j 
during  (t-T-R,t-T]  must  equal  one  of  the  numbers 


0,  Q-(i-k),  Q- ( i-k  )+Q , Q-( i-k  ) + 2Q, . . . . That  is,  the  demand  at  base  j,  n 
must  be  an  element  of  the  set 
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= In:  n = Q-(i-k)  + nQ,  n = 0,1,2,...,  and  n < Iq/Q}  u f0}, 
2 1 

Next,  suppose  I = i < I = k.  Then  demand  at  base  j during 
(t-T-R,t-Tj  must  equal  one  of  the  numbers  k-i,  k-i+Q,  k-i+2Q,...  , 
and,  therefore,  must  be  an  element  of  the  set 


= {n:  n = k-i  + nQ,  n = 0,1,2,...,  and  n <_  I^/Q}. 


We  will  now  find  P(I  = i;  Case  A)  by  considering  separately  the 

2.1  , 

cases  where  I = i > I = k and  l < k.  Thus 


P(  T 2 = i;  Case  A )=  l l l 


iQeJ  k=r+l  neN^ 


i_/Q-[(n-k+Q+r)/Q] 

0 X 

m=0 


P ( T 2 = i | I1=k ; G=m ; D=n;  IQ=i  ) • P( IX=k ; G=m;  D=n ; I0=iQ) 


r+Q  i./Q-[ (n-k+Q+r )/Q] 

III  l 

ip<J  k=i+l  ncN^  m-0 

P( I2=i  l I1=k ; G=m;  D=n;  I =i  ) •P(I1=k;  G=m;  D=n;  I =i  ), 


(Eq.  2) 


1 2 

Observe  that  when  I < I 


P(I2  = i | Ii  = k;  G=m ; D=n ; I0=iQ)  ' 


0,  n i (and  Case  A conditions  hold) 


1,  n c N (and  Case  A conditions  hold); 


• _ 


■HU 
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1 2 

also,  when  I > I 


0,  n i N„ 


P(r=i|li=k;  G=m ; D-n ; lQ=i0)  = 


1,  n e N„ 


(and  Case  A conditions  hold) 


(and  Case  A conditions  hold). 


Furthermore,  due  to  the  independence  of  the  random  variables  1^,  G,  D, 
and  Iq  (see  reference  5), 


P(I  =k;  G=m ; D=n, 


VV  = 


P(I  =k)  • P(G=m)  • P(D=n) 


P,VV' 


But  P( I1=k  ) = 1/Q 


(see  reference  2), 


P(G=m)  = a(-yR,m),  for  m = 0,1,...  (due  to  the  basic  assumptions), 

P(D=n)  = a(AR,n),  for  n = 0,1,...  (since  demand  is  Poisson  distributed 

at  base  j ) , 

and 

P(I  =i.)  = — , for  i„  e J (see  reference  5). 

U U M U 

2 

Hence  P(I  = i;  Case  A)  can  be  found  by  substituting  the  above  probability 
expressions  into  Eq.  2.  Furthermore,  by  substituting  the  resultant 
expression  into  Eq.  1 we  have  obtained  P(B(t)  = b;  Case  A). 

Let  us  now  consider  Case  B and  determine  the  probability  of  observing 
b backorders  at  base  -j  at  time  t given  that  the  total  number  of  units 
demanded  from  the  depot  during  (t-T-R,t-T]  exceeds  the  available  depot 
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inventory;  that  is,  G + [ (D-I  *'+Q+r)/Q  j > Ig/Q.  Observe  that 

B(t ) = max(0,D+U-Q-I2), 

Then 

r+Q  9 9 

P(B(t)  = b;  Case  B)  = £ P(B(t)  = b | Case  B;  I = i)  -P(T  = i;  Case  B) 

i=r+l 

r+Q  [ ( i+b )/Q]  % 2 

= l l P(D=i+b-uQ)  • P(U  = u | Case  B;  I = i) 

i=r+l  u=0 

• P(I2  = i;  Case  B)  for  b ^ 1. 

% o 9 

We  now  find  P(D=i+b-uQ),  P(U  = u|Case  B;  I = i)  and  P(I  = i;  Case  B). 
Since  ft  has  a Poisson  distribution, 

P(ft=i+b-uQ)  = a(AT,i+b-uQ) . 

Rather  than  directly  computing  P(U  = u|Case  B;  I = i),  we  will  find  the 
distribution  of  the  number  of  satisfied  orders  at  time  t that  were  placed 
during  (t-T-R,t-T]  on  the  depot  by  base  j;  that  is,  we  will  determine 
the  conditional  distribution  for  V. 

We  begin  by  evaluating 

P( V=v | D=d ; G=g;  1 1 = k ; IQ  = Qq;  Case  B). 

The  arrival  process  at  base  j is  a Poisson  process;  furthermore,  the 
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process  generating  depot  orders  from  all  other  bases  is  approximately  a 
Poisson  process.  From  base  j's  viewpoint,  system  arrivals  consist  of 
the  demands  placed  at  that  base  and  depot  orders  placed  by  all  other 
bases.  The  sequence  in  which  orders  are  placed  on  the  depot  by  base  j 
and  all  other  bases  determines  the  number  of  satisfied  base  j depot 
orders  at  time  t.  If  D = d and  G = g,  then  every  ordering  of 
the  d+g  arrivals  is  equally  likely  to  occur  since  the  combined  arrival 
process  is  a Poisson  process  (being  the  superposition  of  two  Poisson 
processes--the  base  j arrival  process  and  the  depot  order  process  from 
all  other  bases  ) . 

We  will  compute  P(V=v|D=d;  G=g;  I1=k;  IQ=qQ;  Case  B)  by  considering 
two  cases.  In  the  first  case,  the  last  satisfied  depot  order  is  placed  by 
a base  other  than  j.  Then  if  v orders  are  satisfied  and  I1  = k, 
at  least  vQ-(Q+r-k)  and  not  more  than  vQ-(Q+r-k )+Q-l  demands  must 
occur  at  base  j prior  to  the  arrival  of  the  (IQ/Q-v)-th  order  at  the 
depot  from  a base  other  than  j.  In  the  second  case,  the  last  satisfied 
depot  order  is  placed  by  base  j . 

Combining  the  above  observations  we  see  that 


P(V-v|D=d;  G=g;  I1=k;  IQ=qQ;  Case  B) 


q-1  ( d Y g ) 

T 'svQ-(Q-k+r)+w'Vi-v-r' 
L / 


w=0  ( 

\v0- 


d+g 

vQ- ( Q-k+r )+w+q 


g-(q-v-l) 

g+d- ( vQ- ( Q-k+r )+w+q-v- 1 ) 


(Eq.  3) 


( d )(  g ) 

wQ-( Q-k+r )-y\q-v^  # d-( vQ-(Q-k+r )-l ) 

+ / d+g  "N  g+d- (vQ-(  Q-k+r  )-l+q-v) 

\ vQ-( Q-k+r )+q-v-l/ 
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whenever  v q;  otherwise,  the  probability  is  0.  Also,  we  define 
(^)  = 0 if  z > p,  z < 0,  or  p < 0. 

Now  Lot  us  find 


P(D=d;  G=g ; I1=k;  IQ=qQ ; Case  B). 


This  probability  is  equal  to 


P(D=d;  G=g;  I1=k;  I0=qQ|Case  B;  I2  = i)  . P(Case  B;  I2  = i) 


= P( Case  B ; I2  = i | D=d ; G=g;  IX=k;  I0=qQ)  -P(D=d;  G=g;  IX=k;  IQ=qQ) 


0;  if  [ (d+0+r-k  )/Q]  + g < Iq/Q  = q or  d i wlien  i _> 


1 1 


and  d / N^  when  i < k 


1 • a(AR,d)  • a(yR,g)  • -=■  • tt  ; otherwise. 


Consequent ly , 


P(D=d;  G=g;  I1=k;  Io=qQ|case  B;  I2  = i) 


0;  C(d+Q+r-k)/Q]  + g < q or  d 4 N^  when  i _>  k 

and  d i N^  when  i < k 


va(AR,d)  • a(yR,g)  • i ^ 


Q M P(Case  B;  I2  = i) 


; otherwise 


Now  observe  that  V = v if  and  only  if  U = f (d+0+r-k )/Q]-v  when 


D=d,G=g,I  =k,  and  IQ  = qQ  in  Case  B.  Therefore 
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P(  U=u  | D=d ; C=g ; IJ=k;  IQ=qQ;  Case  B) 

= P(V  = [ (d+Q+r-k)/Q]-u | D=d ; G=g;  I1=k;  IQ=qQ;  Case  B). 

Then  using  Eq.  3 we  may  find  the  conditional  probability  for  the  number  of 
unsatisfied  depot  orders  at  time  t for  orders  placed  by  base  j during 
( t-T-R ,t-T  J . 

Upon  combining  che  previous  results  we  see  that 

P(U=u;  D=d;  G=g;  I^  = k;  IQ=qQ|Case  B) 

= P(U=u|D=d;  G=g;  I^=k;  IQ=qQ;  Case  B)  • P(D=d;  G=g;  I1=k;  IQ=qQ|Case  B), 

and 

r»  / +M 

D/Q  r+Q  ® “> 

P( U=u | Case  B;I=i)=|  l l l 

q=rDy^+l  k=r+l  d=0  g=q-[ (d+Q+r-k )/Q]+l 

P(V  = [ (d+Q+r-  k)/Q]-u | D=d;  G=g;  I1=k;  IQ=qQ;  Case  B) 

• P(D=d ; G=g;  I1=k;  I0=qQ|Case  B ; I2  = i). 

The  final  probabilities  that  we  must  compute  to  complete  the  calculation 

of  P(B(t)  = b;  Case  B)  are  P(I2  = i;  Case  B). 

2 

P(I  = i;  Case  B)  can  be  derived  in  virtually  the  same  way  we  previously 

2 

derived  P(I  = i;  Case  A).  It  is  easy  to  show  that 
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l l l 

k = r+l  ntN^  m= ( i /Q )- [ (n-k+Q+r )/Q ]+l 

~ ■ a(yR,m)  • a(AR,n)  • ~ 

Q M 

r+Q  “> 

l l . I 

k=r+l  ncN^  m= ( iQ/Q )-[ ( n-k+Q+r  )/Q]  + l 
~ • i ■ a(yR,m)  • a( Xp  ,n  ) . 

Thus  we  have  found  P(l3(t)  = b;  Case  A)  and  P(B(t)  = b;  Case  B). 
Since  the  two  cases  form  a partition, 

P(B(t)  = b)  = P(B(t)  = b;  Case  A)  + P(B(t)  = b;  Case  B) 

for  b > 1 and 


P(I2  = i;  Case  B)  = 


l 


i0eJ 


P(B(t)  = 0)  s 1 - l P(B(t)  = b). 
b>l 
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III.  SITUATION  TWO— TWO  BASES  AND  A DEPOT 

In  this  section  we  examine  a system  consisting  of  a depot  and  two 
identical  bases.  The  system's  operation  is  described  in  Section  I, 

Assumptions  1 )— 7 ) stated  in  the  previous  section  are  still  appropriate; 
however,  we  assume  in  this  case  that  there  are  only  two  bases  in  the 
system.  Our  objective,  as  before,  is  to  determine  the  probability 
distribution  for  the  number  of  units  backordered  at  a base  at  an  arbitrary 
point  in  time. 

The  desired  probability  distribution  will  be  developed  by  again  examining 
two  separate  cases.  In  Case  A,  depot  demand  during  (t-T-R,t-T3  is 
assumed  not  to  exceed  IQ  by  time  t-T.  In  Case  B,  depot  demand  exceeds 
I during  (t-T-R,t-T].  Depot  orders  placed  by  a base  may  not  be  satisfied 
by  time  t in  the  latter  case.  To  find  the  probability  distribution  for 
the  number  of  backordered  units  at  a base  at  time  t,  we  examine  the 
arrival  sequences  at  the  bases.  Information  relating  to  the  inventory 
position  at  each  base  and  the  depot  at  time  t-T-R  will  be  combined  with 
knowledge  of  the  arrival  sequence  at  the  bases.  In  the  previous  section, 

G— the  number  of  orders  placed  on  the  depot  by  other  bases— was  Poisson. 

Her-’  it  is  not  and  our  development  parallels  that  of  Section  II  except 
for  this  difference. 

Without  loss  of  generality,  we  will  determine  the  desired  probability 
distribution  for  base  1.  Subscripts  in  this  section  refer  to  a particular 
base.  In  Case  A we  assume  that 

C(D1-I^+Q+r)/Q]  + C(D2-l2+Q+r)/Q]  4 Iq/Q; 
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that  is,  the  total  depot  demand  during  (t-T-R,t-Tj  does  not  exceed  the 
inventory  available  to  meet  that  demand,  1^. 

Since  all  demands  placed  on  the  depot  during  the  period  (t-T-R,t-T] 
are  satisfied  by  time  t in  Case  A,  all  backorders  at  base  1 at  time  t 
are  a result  of  demands  occurring  during  the  interval  (t-T,t].  Thus 

r+Q  -v  9 7 

P(B1(t)  = b;  Case  A)  = £ P(D^  = i+b|Case  A;  1^  = 1)-  Pdf  = i;  Case  A) 

i=r+l 

r+Q  2 

= l a( AT ,i+b ) • P(I  = i;  Case  A),  for  b > 1. 

i=r+l  ~ 

Using  the  same  logic  as  given  earlier  we  see  that 

- i Q+r 

P(I  = i;  Case  A )=  l [ l l l 

iQcJ  k=r+l  neN^  j=r+l  deF 

P(lJ=i|lJ=k;  D2=d;  Dx=n;  I*=j ; I0=iQ) 

•P(lJ=k;  D2=d;  D1=n;  I2  = j;  I0=iQ) 

r+Q  r+Q 

+ l l l II  <Eq.  4) 

iQeJ  k=i+x  neN2  j=r+l  deF 

P(I*=i|lJ=k;  D2=d;  D1=n;  I*= j ; I0=iQ) 

•P(lJ=k;  D2=d;  D^n;  1*= j ; I0=iQ), 

where  and  N?  have  the  same  meanings  as  stated  in  Section  II,  and 

F = {d:  d = 0,1,...;  [(n-k+Q+r)/Q]-Q  + [(d-j+0+r)/Q]-Q  4 iQ}. 


18 


1 2 

But  when  I 1^, 


P<Il  = ilIl=ki  °2=di  Drn;  T21  = j;  W 


1;  n c N ^ (Case  A holds) 


0;  otherwise; 


and  when  I,  > I„, 
x 2 


P(I^i|lJ=k;  D^d;  D1=n;  I*  = j;  W = 


1;  n £ Nj  (Case  A holds) 


0;  otherwise. 


By  independence. 


P(lJ  = k;  D2=d;  Dx=n;  I*=j;  I0=iQ)  = P(lJ=k)  • P(D2  = d)  • PtD^n)  • Pd^j)  • P(I0=i0) 


= — r-  • ^ • a(AR,n)  • a(\R,d). 
Q2  M 


Upon  substituting  these  results  into  Eq.  4 we  have  P(I^  = i;  Case  A). 

In  Case  B,  total  depot  demand  during  (t-T-R,t-T]  exceeds  the 
available  depot  inventory;  that  is. 


f (D1~I1+Q+r)/Q3  + r(D2-I2+Q+r)/Q]  > IQ/0. 


Then 


rtQ  [ ( i+b )/Q] 

P(B  (t)  = b;  Case  B)  = £ £ P(D  =i+b-uQ)  • P(U  = u | Case  B;  H = i) 

’ •'  ^ **■ 


i=r+l  u=0 


Pd^  = i;  Case  B),  b >_  1. 
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% . 2 

Wo  will  now  find  P( D^=i+b-uQ) , P(U^  = u|Case  B;  I,  - i),  and 
r(Tj  - i;  Case  B ) . 

% 

Since  demand  is  Poisson  distributed,  P(D^=i+b-uQ)  = a( XT , i+b-uQ ) . 

. 2 

As  before,  instead  of  determining  P(U^  - u|Case  B;  1^  = i)  directly,  we 
will  first  find  P(  V^v  | f>1=di ; D2=d2;  Il=k;  I2=^’  I0=i0’  CaSe  B)‘  As 
we  saw  earlier,  the  last  satisfied  depot  order  could  have  been  placed  by 
either  base.  Then  conditioning  on  which  base  placed  the  last  satisfied 
depot  order,  we  see  that 


P(V1=v|D1=di;  D2=d2;  I*=k;  I^=j ; I0=i0i  Case  B> 


QyX  /(vQ-(Q+r-k)-l)(vQ  -(Q+r-j  )+w^ 


d^-(vQ-(Q+r-k )-l ) 


“=°  I (< 


Vd2 


v+v)Q-2(Q+r )+k+j 


i+w-1 ^ 


d1+d2~( (v+v)Q-2(Q+r)+k+j+w-l) 


( di  X-  ) 

NWQ-(Q+r-k)+w/>vvQ-(Q+r-  j )-l  Y 

^ / j . j \ 


d2~(vQ-(Q+r-j )-l ) 


f _ dl+d2  ^ d +d  -( (v+v)Q-2(Q+r)tk+j+w-l) 

\(v+v)Q-2(Q+r)+k  + i+w-l  ' 


(Eq.  5) 


where  v = (iQ/Q)-v.  This  result  is  based  on  the  observation  that  each 
sequence  of  total  customer  demands  occurring  at  the  bases — the  super- 
position of  the  two  Poisson  processes — is  equally  likely  to  occur  given 
a fixed  number  of  arrivals  at  each  base. 

Knowing  P(  V ] =v  | D^d^^ ; D2=d2;  I*=k;  1*= j ; I0=i0;  Case  B>  is  equivalent 
to  knowing  the  probability  distribution  for  the  number  ol  uiifiat  inf  i<>.|  depot 
orders  placed  by  base  1 during  (t-T-R,t-Tj. 


Now 


?0 


P(U  = u | Case  B;  I2  = i) 


r+Q  r+Q  i +k-r-l+uQ 

l l l 0 l 


l 


i '‘j  k=r+l  i =r+ 1 d ( =uQ- ( Q+r- i ) d2=iQ- [ ( d L~k+Q+r )/Q ] • 0-(Q+r- j )+ 1 


P(Vi  --  [ (dJ-k+Q+r)/Qj-u|l)i=d1 ; D2  = d?-,  l|=k;  I2  = j;  I0=iQ;  Case 


B) 


P(D1=d1;  U2=d„;  I*=k;  I^= j ; IQ=i0|Case  B;  I2  = i). 


where 


P(D1=d1i  D2=d2’  I"=k’  T2=^’  I0=iolCaSe  B;  T1  = i} 


r a(XR,d.  ) • a(  XR  ,d  ) • • £ — 5 

1 2 02  M PCCase  B;  I?  - i) 


; when  Case  B holds  and 


■ 


d e H when  i >_  k and  d^  e when  i < k 


v.  0 ; otherwise 


(N^  and  have  the  same  interpretation  as  given  in  Section  II), and 

where  Eq . 5 is  used  to  evaluate 


P(VL  = [(d1-k+Q+r)/Q]-u|D1=d1;  D2=d2;  I*=k;  1*= j ; I0=i0;  Case 


B) 


whenever  ' 0;  if  d?  < 0,  this  probability  is  0.  We  will  show  how  to 
2 

compute  P(I^  = i i Case  B). 

following  the  same  reasoning  as  used  previously  for  Case  A,  we  can  show 

that 


1 
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-3  1 r+Q  11 

Pdf  " i*  Case  B)  = l l l l l ^ ^ • a(2XR,n+d) 

IqCJ  k=r+l  ncN^  j=r+l  deH  Q 

r+Q  r+Q  , , 

* l l l l l ^ -i-a(2XR,n+d), 

igtJ  k = i+l  ncN2  j=r+l  d«-H  Q 

where 

H = {d:  d = 0,1,...,  [ (n-k+Q+r )/Q] -Q  + Q(d- j+Q+r )/Q] *Q  > iQ}. 

Upon  combining  the  above  results  we  have  shown  how  to  determine 
P(B^(t)  = b;  Case  B).  As  before,  since  the  two  cases  form  a partition, 

P(B1(t)  = b)  = PCB^t)  = b;  Case  A)  + P(B1(t)  = b;.  Case  B) 

when  b ^ 1 and 

P(B  (t)  = 0)  = 1 - l P(B  (t)  = b). 
b>l 
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IV.  CONCLUDING  COMMENTS 

In  this  paper,  we  have  examined  a two-echelon  inventory  system  in 
which  each  location  follows  a continuous  review  (S,s)  policy.  The 
objective  of  the  paper  was  to  show  how  to  express  the  probability  distri- 
bution for  the  number  of  units  backordered  at  any  base  at  any  arbitrary 
point  in  time.  Two  cases  were  studied.  In  the  first  case,  we  assumed 
there  were  a large  number  of  identical  bases  and  a depot  in  the  system; 
in  the  second  case,  the  system  consisted  of  a depot  and  two  identical 
bases . 

The  results  stated  in  Section  III  can  be  extended  to  the  case 
where  there  are  any  finite  number  of  bases  in  the  system.  We  studied  the 
case  where  there  were  two  bases  only  for  ease  of  exposition.  There  is  no 
theoretical  difficulty  hindering  the  development  of  the  probability  distri- 
butions when  the  number  of  bases  considered  exceeds  two.  Furthermore,  a 
similar  analysis  can  also  be  carried  out  for  the  situation  in  which  the 
A.,  S..  , and  s.  differ  between  the  two  bases.  To  accomplish  this, 

we  drop  assumption  2 and  employ  the  methods  used  in  chapter  5 of  reference  8 
and  the  basic  approach  developed  in  Section  III.  The  probability  distri- 
bution P(B(t)  = b)  is  somewhat  more  tedious  to  calculate  in  the 
i general  case;  however,  there  are  no  theoretical  problems  that  need  to  be 

overcome . 


-V- 1 
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policy.  All  excess  demand  is  assumed  to  be  backordered. 

Our  main  objective  is  to  derive  the  probability  distribution  for  the 
number  of  backordered  units  at  a base  at  an  arbitrary  point  in  time  given 
an  item  follows  a known  (S,s)  policy  at  each  base  and  the  depot.  The 
demand  process  at  each  base  is  assumed  to  be  a stationary  Poisson  process 
The  analysis  is  carried  out  for  two  cases.  In  the  first  case,  we  assume 
the  system  consists  of  a large  number  of  bases;  in  the  second  case,  we 
assume  there  are  two  bases  in  the  system.  To  simplify  the  discussion,  we 
assume  in  both  cases  that  all  bases  follow  the  same  (S,s)  policy. 
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